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Abstract 

We show how to use operators in the description of exchanging processes often taking 
place in (complex) classical systems. In particular, we propose a set of rules giving rise to 
an hamiltonian operator for such a system S, which can be used to deduce the dynamics 
of S 



I Introduction and motivations 



In a series of recent papers we have used an operatorial approach in the description of classical 
systems, with few or with many degrees of freedom, pQ-[8]. In particular, we have shown 
how canonical commutation and anticommutation relations (CCR and CAR respectively) can 
be used in the analysis of simplified stock markets, as well as in the description of simpler 
dynamical systems, like those arising from love affairs. We have also adopted the same general 
settings in the analysis of migration processes and of population dynamics. 

The main ingredient in our approach is the hamiltonian operator H of the system S we 
are interested in, which is used to deduce the time evolution of S, see below. This paper is 
devoted to discuss a minimal set of rules which should be adopted to write down H. Some 
examples of hamiltonians found this way will be discussed. However, the dynamical content of 
these hamiltonians will not be considered here, since it was already discussed elsewhere, p]-[E]. 

The paper is organized as follows: in the rest of this section we review few known fact on 
CCR. We will not discuss here CAR since they will only play a minor role in Section III. 

In Section II we propose our set of rules useful to determine the analytic expression of the 
hamiltonian of a system S. 

In Section III we discuss few examples, while Section IV contains our conclusions. 

The reason why the operator H assumes a crucial role in our approach is because the 
dynamical behavior of S is here assumed to be given by the Heisenberg equation of motion: 

let % be an Hilbert space and B(Ji) the set of all the bounded operators on Let S be 
our physical system and 21 the set of all the operators useful for a complete description of S, 
which includes the observables of S. The time evolution of S is assumed to be given by the 
self-adjoint hamiltonian H = H ' of S, which in standard quantum mechanics represents the 
energy of S. In the Heisenberg picture the time evolution of an observable X e 21 is expressed 



where [A, B] := AB — BA is the commutator between A and B. The time evolution defined 
in this way is usually a one parameter group of automorphisms of B(7i). It might seem that 
adopting the Heisenberg picture in the description of classical systems may appear unappro- 
priate. However, as discussed in our previous literature as well as in many other papers on 




i[H,X{t)], 
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similar subjects, see for instance [9]- [2], this approach is justified a posteriori since, at least for 
simple systems, it produces exactly that time evolution which one expects to find. We should 
also mention that the uncertainty principle arising from the non abelianity of the operators 
involved in the description of S, does not appear in our approach, since all the observables of 21 
do commute. Other authors, on the other hand, because of what they are interested in, consider 
such an uncertainty a richness and not a problem of a quantum view to complex systems, pj)] . 

In our approach a special role is played by the so called CCR: we say that a set of operators 
{a;, a], I — 1, 2, . . . , L) satisfy the CCR if the following hold: 

[ai, a~l] = 5i n l, [ai, a n ] = [a], a*] = 0, (1.3) 

for all l,n — 1, 2, . . . , L. Here 1 is the identity operator on "H. These operators, which are widely 
analyzed in any textbook in quantum mechanics, see [16] for instance, are those which are used 
to describe L different modes of bosons. From these operators we can construct hi = ajai and 
N = Yli=i ™i which are both self-adjoint. In particular hi is the number operator for the 1-th 
mode, while N is the number operator of S. 

The Hilbert space of our system is constructed as follows: we introduce the vacuum of the 
theory, that is a vector <£>o,o,...,o which is annihilated by all the operators af. a/<£>o,o,...,o = for 
all I — 1, 2, . . . , L. Then we act on y2o,o,....o with the operators a\ and their powers: 

<p ni , n2 ,..., nL := - 1 (al)" 1 (4r " " " (4r^o,o,..,o, (1-4) 

Vni! n 2 \ . . . n L \ 

ni = 0, 1,2, . . . for all I. These vectors form an orthonormal set and are eigenstates of both 
hi and N: hiip nun2 ,... t n L = ni<p nit n 2 ,... t n L and N<p nun2 ^ nL = N(p nitn2t _ >nL , where N = Ya=i u i- 
Moreover using the CCR we deduce that 

hi (atyni,n2,...,n L ) = ( n l ~ l)(°Wn 1 ,n 2r ..,nJ (I- 5 ) 

and 

hi {AVniw,...^ = fa + 1 )( a j ( Pn 1 ,n 2 ,...,n L ), (1-6) 

for all I. For these reasons the following interpretation is given: if the L different modes of 
bosons of S are described by the vector ^p ni ,n 2 ,...,n L i this implies that ri\ bosons are in the first 
mode, 712 in the second mode, and so on. The operator hi acts on f ni ,n 2 ,...,n L and returns 
ni, which is exactly the number of bosons in the 1-th mode. The operator iV counts the total 
number of bosons. Moreover, the operator ai destroys a boson in the 1-th mode, while a\ creates 
a boson in the same mode. This is why ai and a\ are usually called the annihilation and creation 
operators. 
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The Hilbert space T-L is obtained by taking the closure of the linear span of all these vectors. 
A similar construction can be repeated starting with CAR, but we will not consider this 
possibility here since it is not essential for the general analysis we will discuss in this paper. 

An operator Z £ 21 is a constant of motion if it commutes with H . Indeed in this case 
equation (11.21) implies that Z(t) = 0, so that Z(t) = Z for all t. 

The vector (p nij n 2 ,.~,nL * n P-^P defines a vector (or number) state over the algebra 21 as 

>> (1-7) 

where ( , ) is the scalar product in T-L. As we have discussed in [1]- [8], these states are used to 
project from quantum to classical dynamics and to fix the initial conditions of the system. 



II The rules 

As already discussed in the Introduction, the main interest in this paper is to produce a sort 
of recipe which has to be used to write down the hamiltonian H of the classical system S we 
are interested in. To simplify our analysis, let us first suppose that S consists of two main 
interacting parts, S\ and S2 , the actors of the game, whose union reproduces S and which have 
no intersection: S = Si U S2 and S\ D S2 = 0. Suppose now that S\ and S2 can exchange 
something, OJt, which can only take integer values^. A typical example of this situation is in 
stock markets, where two traders Si and S2 , exchange money (and shares). Other examples 
are discussed in [6] and [7j, where what is exchanged is mutual affection (in other words, love!) 
between the two lovers. In [8] we have two populations in different regions of a two-dimensional 
lattice, and they exchange people, i.e. there is people moving from one region to the other. 

Let us now introduce two annihilation operators, a\ and 02, related respectively to DJli 
and 30^2 5 and their conjugate creation operators a\ and a\. Here OJTi is that part of QJt which 
belongs to S\: the money of the first trader, or the number of shares in his portfolio, or jet 
the amount of love that Bob, the first lover, experiences for Alice, the second one, and so 
on. As in the Introduction, these operators obey the following CCR: [a$, c^] = [a|,at] = 0, 
[aj, at] = 5ij 1. Calling <£> ,o the vacuum of a±, a 2 , that is that vector of T-L annihilated by a\ 

1 \ n \ \ n 2 i >i ... 

and 02, 01^0,0 = a 2V?o,o = 0, the vector <f ni>n2 '■= v / Wl i TO2 i a \ a 2 ¥0,0 describes a situation m 
which the value of VJl\ is n\ and that of 9Jt2 is n2\ indeed, calling hj := OjOj the related number 
operators, we know that hjLp ni ^ n2 = rij <p nun2 , j = 1,2. In all the papers written so far, p]-[E], 

^^We could relax this assumption by assuming that the values of 9JI are discrete rather than integer. 
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the eigenvalues of fij, rij, are considered directly related to the value of 9R j. Now we are ready 
to state the first rule of our construction: 



Rule 1:—The exchange ofDJt between S\ and S 2 is modeled adding to the hamiltonian of S 
a term a\ 02 + a\ 0,1 ■ If, for some reason, the model should be non-linear, then this contribution 
must be replaced by a\ a 2 + a\ af 1 , M > 1 being a measure of the non linearity. 

The motivation of this rule is given by the action of a\ a 2 + a\ a>i on the vector </ 3 m,n 2 : 



where the normalization constants are missing since they are not interesting for us, here. As 
we can see, what we get is a combination of two vectors: the first one, V9 ni+l n2 _ 1 , shows that 
the value of DJti is increased by one unit while, simultaneously, DJl 2 decreases by a unit. In 
the second contribution, (p ni _ 1>n2+1 , the opposite happens. In both cases, what it is going on 
is that Si and S 2 are exchanging one unit of Wl. Analogously, acting with a\ a 2 + a\ af 1 
on ifini,n 2 would produce a combination of vectors f ni+ M,n 2 -i an d ^ m -M,n 2 +i) which is useful 
to introduce a possible asymmetry between S± and S 2 , or, from a dynamical point of view, a 
non-linearity in the dynamics of S, [6]. 

One may argue why not to add simply a\ a 2 in the hamiltonian. The reason is the following: 
if we don't consider both a\a 2 and a\ai, the final hamiltonian would not be self-adjoint, and 
this would create a lot of difficulties in finding a reversible time evolution: for instance, if 
H 7^ then, among other problems, the norm of e lHt Xe~ lHt is different from that of X, so 
that the probabilistic interpretation of the wave-function typical of quantum mechanics would 
be lost. This is obviously related to the decay effects which we don't want to discuss here. We 
will briefly come back on this aspect in Section IV. 

Let us now go the the second rule of our construction: 

Rule 2:~ The hamiltonian H for S must contain a term, Hq, such that, in absence of 
interaction between S\ and S 2 , their related number operators, hi and h 2 , stay constant in 
time. 

This is quite a natural assumption: if S\ and S 2 do not interact, there is no reason for 
them to modify their situation, and in particular there is no reason (and no possibility!) for 
exchanging units of 9JT. To be concrete, this means that, if at t = S is described by the 
state (p ni ,n 2 , and if no interaction is contained in the hamiltonian, H = Hq, then at t > the 
system is still described by (p ni ,n 2 (but, at most, for an overall phase). It should be stressed 
that, however, this does not imply that in this case there is no dynamics at all! What we are 
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claiming is that hj(t) = hj(0), but this does not imply that, for instance we also necessarily 
have aj(t) = %(0). On the contrary, in many examples this is not so, see [T] for such an 
example. 

There is still another rule which is quite useful in the determination of H. For that it may 
be convenient to recall the notion of closed system: a system S is called closed if it has no 
interaction with the environment 1Z. 

Rule 3:— If S is a closed system, the hamiltonian H of S must commute with those global 
number-like operators related to the observables which are not exchanged between S and 1Z. 

The motivation is, again, rather natural: as we have seen in the Introduction, all the 
observables which commute with the hamiltonian are integrals of motion, so that they do not 
change with time. This is exactly what is expected to the global quantities of the system S, 
since they are not moving outside S. A simple example of this situation is provided by the 
total number of shares of a certain type in a closed market where the shares are not created 
or destroyed: if at t = this number is N = ri\ + n 2 , where rij is the number of shares of that 
kind which belong to Sj, then N does not change with time, even if the number of shares in 
each trader's portfolio does change, in general. In this case, Rule 3 reads \H,hi + h 2 ] = 0, 
while [H, hi] ^ and [H, h 2 ] ^ 0, in general. 

Ill Examples 

In this section we will show how the rules described so far can be explicitly used in the analysis 
of some classical systems, and which kind of hamiltonian are deduced. 

III.l First example: love affair 

The first model we have in mind consists of a couple of lovers, Bob and Alice, which mutually 
interact exhibiting a certain interest for each other. Of course, there are several degrees of 
possible interest, and to a given Bob's interest for Alice (LoA, level of attraction) there corre- 
sponds a related reaction (i.e., a different LoA) of Alice for Bob. In our previous decomposition 
of S in Si and S 2 , here Bob plays the role of Si, while Alice that of S 2 , and 971 is the mutual 
affection between the two. The bosonic operators associated to Bob are ai, a[ and hi = a\ai, 
while those associated to Alice are a 2 , a\ and h 2 = a\a 2 . The (integer) eigenvalue rii of hi 
measures the value of the LoA that Bob experiences for Alice: the higher the value of ni the 
more Bob desires Alice. For instance, if n± = 0, Bob just does not care about Alice. We use 
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n 2 , the eigenvalue of h 2 , to measure the attraction of Alice for Bob. The law of attraction we 
have in mind states that, if n\ increases, then n 2 decreases and viceversa. This suggests to use 
the following self-adjoint operator to describe the interaction between Alice and Bob: 



H = \(a¥4 + a 2 a\ M y (3.1) 



where M describes a sort of relative behavior, [6]. This choice is written following Rule 1 of 
the previous section and it trivially satisfies Rule 2: if A = there is no dynamics at all since 
H = and, as a consequence, [H, hi] = [H, h 2 ] = 0. Concerning Rule 3, it is an easy exercise 
to check that I(t) := hi(t) + is a constant of motion: I(t) = 1(0) = h\(0) + M h 2 (0), 

for all t 6 K, since [H, I] = 0. Therefore, during the time evolution, a certain global attraction 
is preserved and it can only be exchanged between Alice and Bob: notice that this reproduces 
our original point of view on the love relation between Alice and Bob: the more Bob falls in 
love with Alice, the less Alice cares about Bob! If M is fixed to be one then Bob and Alice 
react in the same way and the model becomes exactly solvable, [6]. 

In [6] we have also considered a love affair involving, other than Alice and Bob, a third 
actress, Carla, also having a relation with Bob. Our assumptions are the following: (1) Bob 
can interact with both Alice and Carla, but Alice (respectively, Carla) does not suspect of 
Carla's (respectively, Alice's) role in Bob's life; (2) if Bob's LoA for Alice increases then Alice's 
LoA for Bob decreases and viceversa; (3) analogously, if Bob's LoA for Carla increases then 
Carla's LoA for Bob decreases and viceversa; (4) if Bob's LoA for Alice increases then his LoA 
for Carla decreases (not necessarily by the same amount) and viceversa. 

Introducing now the operators a 3 , a\ and n 3 = a\a 3 for Carla, and splitting the operators 
related to Bob in two (i.e. a±2 and 013 to describe the interaction between Bob and, respectively, 
Alice and Carla) the hamiltonian which describes all these effects is the following: 

H = Ai 2 (a{ 2 a 2 + ai2 4) + A i3 ( a i3 a 3 + 013 oJJ + Ai (a\ 2 a 13 + a 12 oj 3 J , (3.2) 



for some real values of A's. It can be easily seen that the first contribution, A12 ya\ 2 a 2 + ai2 4 

describes the mechanism (2) above, while A13 (^a\ 3 a% + 013 a\ J is related to point (3). Point (4) 

is implemented by Ai (a\ 2 ais + ai2 a{ 3 j . These three contributions all trivially satisfy Rule 
1 and, again, Rule 2 is also verified: no interaction means that all the A's are zero, so that 
H reduces to the zero operator, and all the observables stay constant in time. Let us now 
introduce h\ 2 = a\ 2 a\ 2 , describing Bob's LoA for Alice, hi 3 = a\ 3 a\ 3 , describing Bob's LoA 
for Carla, h 2 = a\a 2 , describing Alice's LoA for Bob and h 3 = a\a 3 , describing Carla's LoA 
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for Bob. If we define J := h\% + his + ^2 + ^3, which represents the global level of LoA of 
the triangle, this is a conserved quantity: J(t) = J(0), since [H, J] = 0: no exhange with the 
environment is possible, here! It is also possible to check that [H, h\i + 7T. 13 ] 7^ 0, so that the 
total Bob's LoA is not conserved during the time evolution. 

More details, as the equations of motion arising from these hamiltonians, their solutions 
and more extensions can be found in [6] and [7]. 

III. 2 Second example: competition between species and migration 

In this example we consider a two-dimensional region TZ in which two populations S\ and S2 are 
distributed. In [8] we have considered these species as predators and preys, or as two migrant 
populations, moving from one part of TZ to another. Following the above rules we can construct 
the hamiltonian of the full system S. However, for reasons discussed in [8], it is convenient to 
use here annihilation and creation operators satisfying CAR rather than CCR. This choice is 
motivated by a first technical and a second more substantial reason: the technical reason is that 
we get finite dimensional Hilbert space for S, while the more substantial reason is that, using 
CAR, we automatically incorporate an upper bound for the densities of the two populations, 
which is a natural requirement for our biological interpretation. 

The starting point is the (e.g., rectangular or square) region TZ, which we divide in iV cells, 
labeled by a = 1, 2, . . . , N. In each cell a the two populations, whose related operators are a a , 
a a and ha = o) a a a for what concerns «Si, and b a , b a and ha = b^ a b a for iS 2 , are described by 

H a = H° a + \ a Hl, H° a = u a a ala a + u b a blb a , H J a = a\b a + b\a a . (3.3) 

It is natural to interpret the mean values of the operators h^a and nh} as local density operators 
of the two populations in the cell a: if the mean value of, say, n£\ in the state of the system is 
equal to one, this means that the density of S\ in the cell a is very high. Notice that H a = iff, 
since all the parameters, which in general are assumed to be cell-depending (to allow for the 
description of an anisotropic situation), are real and positive numbers. The CAR are 

{a a ,al} = {b a ,bl} = 6 atfi l, {<4&Jj} = 0. (3.4) 

Of course, the full hamiltonian H must consist of a sum of all the different H a plus another 
contribution, Hdiff, responsible for the diffusion of the populations all around the lattice. A 
natural choice for H^ff, in view of the above rules, is the following: 

Hdiff = ^2p a ,p 1 7a (a> a a\ + a>pa a ) + lb (b a b ] p + bpb a J j , (3.5) 
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where also 7 a , 7& and the p a ^ are real quantities. In particular, p a ^ can only be or 1 depending 
on the possibility of the populations to move from cell a to cell (3 or vice-versa. For this reason 
they are considered as diffusion coefficients. Notice that a similar role is also played by j a and 
7 fe . H = ^2 a H a + H di ff obeys the three rules of Section II. Indeed, if there is no interaction 
between Si and S2, and between members of the same species localized in different cells of 7Z, 
it is easy to check that the densities of Si and S2 stay constant in all the cells: ha {t) = (0) 
and ha\t) = ha\o), for all a. Hence Rule 2 holds true. Rule 3 is also satisfied, since Si and 
1S2 cannot move outside 7Z: it is again possible to find an operator, related to the total number 
of members of S distributed all along 7Z, which commutes with H, so that this global density 
stays constant in time. Concerning Rule 1, we see that this is applied several times in the 
definition of H. For instance we have the contribution a) a b a + b a a a , which shows how Rule 1 
is applied in the interaction between Si and S 2 in the cell a, but we also have a a a^ + apa a , 
which is again Rule 1, but applied to Si in different cells. And so on. Again, we refer to [8] 
for the analysis of the equations of motion arising from this hamiltonian. 

III. 3 Last example: stock market 

In recent years we have proposed several hamiltonians describing simplified stock markets, 
[T]-j4]. The one we discuss here, the most efficient proposal, so far, was first introduced in [4]. 

Let us consider iV different traders T\, T2, . . ., tn, exchanging L different kind of shares 
o~i, cr 2 , . . ., ai. Each trader has a starting amount of cash, which is used during the trading 
procedure: the cash of the trader who sells a share increases while the cash of the trader who 
buys that share consequently decreases. The absolute value of these variations is the price of 
the share at the time in which the transaction takes place. It is clear that the above-mentioned 
division of S in just two plus one components, «Si, S2 and 97t, must be extended here, while 
the main ideas are unchanged. It is convenient to introduce a set of bosonic operators which 
are listed, together with their economical meaning, in the following table. We are adopting 
here latin indexes to label the traders and greek indexes for the shares: j = 1,2, . . . , N and 
a = 1,2, ...,L. 
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the operator and.. 


...its economical meaning 






annihilates a share a a in the portfolio of re- 






creates a share a a in the portfolio of re- 






counts the number of share a a in the portfolio of Tj 






annihilates a monetary unit in the portfolio of re- 




4 


creates a monetary unit in the portfolio of re- 




^3 = C j C j 


counts the number of monetary units in the portfolio of Tj 




Pa 


lowers the price of the share a a of one unit of cash 




Pa 


increases the price of the share a a of one unit of cash 




Pa PaPoi 


gives the value of the share a a 



Table 1.- List of operators and of their economical meaning. 



These operators are bosonic in the sense that they satisfy the following commutation rules 

\ C V 4] = 1 $3,k> \Pa,Pp] = 1<W [ a j,OL,a\^\ = l<^,fc<W, (3-6) 

while all the other commutators are zero. We assume that the hamiltonian of the market, H, 
can be written as H = H + H pr i ces , where 

H = H + XH t , with 
< Ho = Y.j,a u j,a + 12 j u j k j (3.7) 
Hi = V,. M! pg> ^.acf-c}^ + ft.c.) . 

Here h.c. stands for hermitian conjugate, cf a and cj " are defined as in [2], and ojj, a , ujj and p-^ 
are positive real numbers. In particular these last coefficients assume different values depending 
on the possibility of r, to interact with Tj and exchanging a share o~ a : for instance = if 
there is no way for r 2 and r 5 to exchange a share o\. It is natural to put p\°^ = and p-^ = pj°p. 

Going back to ( 13. 7ft . we observe that if obeys Rule 2 of Section 2, since, if there is no 
interaction between the traders, then A = and, as a consequence, H = Hq: [Hq, hj, a } = 0, for 
all j and a. As for Hi, this is written obeying Rule 1: the action of a single contribution of Hi, 
a l,a a j,a c i' ac j a i on a vec tor number which extends those introduced in Section II, (p^ rij a }-{k 3 };{p a }i 
is proportional to another vector ip{ n >. };{k'.};{P^} with just 4 different quantum numbers. In 
particular rij >a , n ija , kj and ki are replaced respectively by rij >a — 1, n^ a + 1, kj + P a and fcj — P a 
(if this is larger or equal than zero, otherwise the vector is annihilated). This means that Tj is 
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selling a share o~ a to Ti and earning money from this operation. For this reason it is convenient 
to introduce the following selling and buying operators: 

x j,a '■= a j,a c ] ) x j, a := a ],a C j " (3-8) 

With these definitions and using the properties of the coefficients pfj we can rewrite Hi as 

H i = 2 Yl PiJ x U X J^ H = Yl u i' a h i> a + Y u i~ k i + 2X Y1 Pi°3 )x U x h« ( 3 - 9 ) 

i,3,a j,a j i,j,a 

The role of H pr i Ce in [I] was to fix the time evolution of the operators P a , a = 1, 2, . . . , L. The 
hamiltonian H corresponds to a closed market where the money and the total number of shares 
of each type are conserved. Indeed, calling N a := X^=i ^i,a an d K '■= J2iLi h we see that, for 
all a, [H, N a ] = [H, K\ = 0. Hence N a and K are integrals of motion, as expected: Rule 3 is 
satisfied. We refer to jl] for the analysis of the time evolution of the portfolio operator of the 
trader r h n,(t) = YlLi ^(t) m, a (t) + fc,(t). 



IV Further considerations and conclusions 

The same rules have already been adopted for systems which are not closed, i.e. for those 
systems which exchange something with the environment. This is discussed, for instance, in 
[7]: again, the main idea is that we can use creation and annihilation operators also in the 
description of the reservoir, and in modeling an exchange between the system and the reservoir. 
This exchange is described adding in the hamiltonian a contribution obeying Rule 1, while 
Rule 2 has to be intended here in the following way: if S does not interact with the reservoir, 
no decay is allowed. Rule 3 is recovered for some global quantity which mixes the degrees of 
freedom of the reservoir and of the system. The conclusion is, therefore, that our rules can be 
used also in more general, and sometimes more useful, contexts. 
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